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ABSTRACT 

The backreaction of the Lorentz force on the a-effect is studied in the limit of small magnetic 
and fluid Reynolds numbers, using the first order smoothing approximation (FOSA) to solve 
both the induction and momentum equations. Both steady and time dependent forcings are 
considered. In the low Reynolds number limit, the velocity and magnetic fields can be ex- 
pressed explicitly in terms of the forcing function. The nonlinear a-effect is then shown to be 
expressible in several equivalent forms in agreement with formalisms that are used in various 
closure schemes. On the one hand, one can express a completely in terms of the helical prop- 
erties of the velocity field as in traditional FOSA, or, alternatively, as the sum of two terms, 
a so-called kinetic a-effect and an oppositely signed term proportional to the helical part of 
the small scale magnetic field. These results hold for both steady and time dependent forcing 
at arbitrary strength of the mean field. In addition, the r-approximation is considered in the 
limit of small fluid and magnetic Reynolds numbers. In this limit, the r closure term is absent 
and the viscous and resistive terms must be fully included. The underlying equations are then 
identical to those used under FOSA, but they reveal interesting differences between the steady 
and time dependent forcing. For steady forcing, the correlation between the forcing function 
and the small-scale magnetic field turns out to contribute in a crucial manner to determine the 
net a-effect. However for delta-correlated time-dependent forcing, this force-field correlation 
vanishes, enabling one to write a exactly as the sum of kinetic and magnetic a-effects, similar 
to what one obtains also in the large Reynolds number regime in the r-approximation closure 
hypothesis. In the limit of strong imposed fields, Bo, we find a oc Bq for delta-correlated 
forcing, in contrast to the well-known a oc Bq 3 behaviour for the case of a steady forcing. 
The analysis presented here is also shown to be in agreement with numerical simulations of 
steady as well as random helical flows. 
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1 INTRODUCTION 

Turbulent mean field dynamos are thought to be at the heart of mag- 
netic field generation and maintenance in most astrophysical bod- 
ies, like the sun or the galaxy. A particularly important driver of the 
mean field dynamo (MFD) is the a-effect which, in the kinematic 
regime, depends only on the helical properties of the turbulence. It 
is crucial to understand how the a-effect gets modified due to the 
backreaction of the generated mean and fluctuating fields. 

Using closure schemes or the quasi-linear approximation it 
has been argued that, due to Lorentz forces, the a-effect gets 
"renormalized" by the addition of a term proportional to the cur- 
rent helicity of the generated small scale magnetic fields (Pouquet 
et al. 1976; Gruzinov & Diamond 1994; Kleeorin et al. 1995; Sub- 
ramanian 1999; Kleeorin & Rogachevskii 1999; Blackman & Field 
2000; Radler et al. 2003; Brandenburg & Subramanian 2005a). The 
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presence of such an additional term is uncontroversial if a helical 
small scale magnetic field is present even in the absence of a mean 
field. However, it has been argued that, in the absence of such a pre- 
existing small scale magnetic field, the a-effect can be expressed 
exclusively in terms of the velocity field, albeit one which is a so- 
lution of the full momentum equation including the Lorentz force 
(Proctor 2003; Radler & Rheinhardt 2007). In the latter case, it is 
not obvious that the helicity of the small scale magnetic field plays 
any explicit role in the backreaction to a. It is important to clarify 
this issue, as it will decide how one should understand the satura- 
tion of turbulent dynamos, as well as the possibility of catastrophic 
quenching of the a-effect and ways to alleviate such quenching. 
Here and below, "catastrophic" means that a is quenched down to 
values on the order of the inverse magnetic Reynolds number. 

In order to clarify these conflicting views, we examine here 
an exactly solvable model of the nonlinear backreaction to the a- 
effect, where we assume small magnetic and fluid Reynolds num- 
bers. Obviously, this approach does not allow us to address the 
question of catastrophic quenching of astrophysical dynamos di- 
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rectly, but it allows us to make novel and unambiguous statements 
that help clarifying the nature of magnetic saturation. We will show 
that, at least in this simple context, both the above viewpoints are 
consistent, if interpreted properly. 



2 MEAN FIELD ELECTRODYNAMICS 

In mean field electrodynamics (Krause & Radler 1980; Moffatt 
1978), any field F is split into a mean field F and a 'fluctuating' 
small scale field /, such that F = F + f. The fluctuating veloc- 
ity (or magnetic) field is assumed to possess a correlation length I 
small compared to the length scale L of the variation of the mean 
field. The magnetic field obeys the induction equation, 

BB 

— =j?V 2 B + V x (U x B), V-B = 0, (1) 

where U represents the fluid velocity, rj = (/i o-) _1 is the mag- 
netic diffusivity (assumed constant), a is the electric conductivity, 
and fio is the vacuum permeability. Averaging Eq. (1), we obtain 
the standard mean-field dynamo equation 

— =?7V 2 B + V x (U x B + S), V-B = 0. (2) 

This averaged equation now has a new term, the mean electromo- 
tive force (emf) £ — u x b, which crucially depends on the statis- 
tical properties of the u and b fields. The central closure problem 
in mean field theory is to find an expression for the correlator £ in 
terms of the mean fields. 

To find an expression for £ , one needs the evolution equations 
for both the fluctuating magnetic field b and the fluctuating velocity 
field u. The first follows from subtracting Eq. (2) from Eq. (1), 

^ =r? V 2 6 + V x (Ux b + u xB)+G, V • b = 0. (3) 
at 

Here G = V x (u x b)' with (u x b)' = u x b — u x b. In what 
follows, we will set the mean field velocity to zero, i.e. U = and 
focus solely on the effect of the fluctuating velocity. 

The evolution equation for u can be derived in a similar man- 
ner by subtracting the averaged momentum equation from the full 
momentum equation. We assume the flow to be incompressible 
with V • u = 0. We get 

% = ^( P +±B-b)+vV 2 u 

at P \ "° / 

+ — — \{B ■ V)6 + (b ■ V)Sl + f + T. (4) 
pap 1 1 

Here p is the mass density, p is the perturbed fluid pressure, v is the 
kinematic viscosity taken to be constant, / is the fluctuating force, 
and 

T = — (u • Vu)' - — [(&■ V6)' - ±V(6 2 )'l (5) 
Pop 

contains the second order terms in u and b. Here, primed quantities 
indicate deviations from the mean, i.e. X' — X — X. We will also 
redefine b/^/Jbop — > b, by setting pop = 1, so that the magnetic 
field is measured in velocity units. 

In order to find £ under the influence of the Lorentz force 
one has to solve Eqs (3) and (4) simultaneously and compute 
u x b. In general this is a difficult problem and one has to take 
recourse to closure approximations or numerical simulations. To 
make progress we assume here R m = ul/r) <C 1 and Re = 
ul/v <§C 1; that is both the magnetic and fluid Reynolds numbers 



are small compared to unity. In this case there is no small scale 
dynamo action and so the small scale magnetic field is solely due 
to shredding the large scale magnetic field. Here u and b (see be- 
low) are typical strengths of the fluctuating velocity and magnetic 
fields respectively. In the low magnetic Reynolds number limit the 
ratio of the first nonlinear term in G to the resistive term in Eq. (3) 
is ~ (ub/l)/(r]b/l 2 ) ~ i? m < 1. So this part of G can be ne- 
glected compared to the resistive term. (Note that the second term 
in G vanishes automatically when taking the averages to evaluate 
the mean emf.) Neglecting the nonlinear term, the generation rate 
of b is ~ uB/l, while its destruction rate is ~ rjb/l 2 . Equating 
these two rates, this also implies that b ~ R m B and the fluctu- 
ation field is only a small perturbation to mean fields. Similarly 
the ratio of the nonlinear advection term to the viscous term in 
Eq. (4), is ~ (u 2 /l)/(vu/l 2 ) ~ Re < 1 and the ratio of the 
parts of the Lorentz force nonlinear in b to that linear in b is 
~ (b 2 /l)/(b~B/l) ~ fl m < 1. So T can also be neglected in 
Eq. (4). In this limit, one can therefore apply the well known first 
order smoothing approximation (FOSA). It is sometimes also re- 
ferred to as the second order correlation approximation (or SOCA; 
see, e.g., Krause & Radler 1980). This approximation consists of 
neglecting the nonlinear terms G and T, to solve both the induc- 
tion and momentum equation. Since FOSA is applied to the mo- 
mentum equation as well, we will refer to this as "double FOSA". 
In order to make the problem analytically tractable, we will take 
B = Bo = const. This also allows us to isolate the a-effect in a 
straightforward fashion. In the next section we begin by consider- 
ing for simplicity the case of steady forcing. It is then possible to 
also neglect the time derivatives in Eqs (3) and (4). We return to 
consider time dependent forcing in detail in Section 4. 



3 COMPUTING £ FOR STEADY FORCING 

Under the assumptions highlighted above, one can solve directly 
for u and b in terms of the forcing function /. This in turn allows 
the calculation of the mean emf in four ways. 

A. We use the induction equation to solve for b in terms of u. 
Then one can write the emf completely in terms of the velocity 
field, as in normal FOSA and then substitute for u in terms of /. 

B. We compute £ = u x b directly. 

C. We use the momentum equation to solve for u in terms of b 
and the forcing function /, and then substitute for b in terms of /. 

D. Compute £ from the d£/dt = relation, as in r- 
approximation closures. 

We will show that all four methods give the same answer for the 
mean emf in terms of the forcing function /. The first Method A 
gives the traditional FOSA result for the a-effect being dependent 
on the helical properties of the velocity field, while Method C can 
be interpreted to reflect the idea of a renormalized a due to the 
helicity of small scale magnetic fields. But we show that the final 
answer in terms of the forcing is identical. 

Before going into the various methods as highlighted above, 
we solve for u and b in terms of the forcing function /. In the low 
conductivity limit, neglecting the time variation of b in Eq. (3) we 
have, 

-r?V 2 b = Bo ■ Vu. (6) 
Similarly, in the limit of low Re and R m , Eq. (4) becomes, 
-vV 2 u = B Vb + f - Vp off , (7) 
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where p e ff combines the hydrodynamic and the magnetic pressure. 
Using the incompressibility condition, one can eliminate p e g . We 
will solve these equations in Fourier space. Throughout this paper 
we will be using the convention 



which satisfies the inverse relation 

u(x) = J ii(k)e ikx dk. 

In Fourier space, Eqs (6) and (7) become 
r]k 2 bi(k) = (ifc • B )ui(k), 
vk 2 Ui{k) = (ifc ■ Bo) bi(fc) + fi(k), 



(8) 

(9) 

(10) 
(11) 



where we have chosen the forcing to be divergence free, with 
ifc • / = 0. We can therefore solve the above two equations si- 
multaneously to express u and b completely in terms of /, 



Ui(fc) 
bi(k) ■ 



Mk) 



vk 2 + (Bo • k) 2 /r)k 2 ' 

fi{k) ifc -Bp 

vk 2 + (B ■ k) 2 /rjk 2 r)k 2 ' 



(12) 



(13) 



We can use these solutions to calculate £. For getting an explicit 
expression, we also need the equal time force correlation function. 
For isotropic and homogeneous forcing, this is given by 



(14) 



fj(P,t) fk(q,t) = 8 s (p + q)F jk (q). 
Here, Fj k is the force spectrum tensor which is given by 

M ' - jk 4nk 2 3fem 87rfc 4 ' 

where Pj k = 8j k — kjkk/k 2 is the projection operator, and <£>(fc) 
and x(k) are spectra characterizing the mean squared value and the 
helicity of the forcing function, normalized such that 



(15) 



$(fc) dfc = \f = \A\. 



f 

Jo 



X (k) dk = f ■ (V x /) = H { . 



(16) 
(17) 



The mean emf can be written as 



£i(x) = e ijk Uj(x)b k (x) = J £i{k) e x dfc, (18) 
where the Fourier transform £ is given by 

£i(k) = e l:jk J Uj(k- q)b k (q) dq. (19) 

We now turn to the calculation of the nonlinear mean emf and the 
resulting nonlinear a-effect in the four different methods outlined 
above. 



At this stage one can put the emf completely in terms of the velocity 
field and recover the usual FOSA expression that in the low con- 
ductivity and isotropic limit, the a-effect is related to the helicity of 
the velocity potential (Krause & Radler 1980; Radler & Rheinhardt 
2007). This can be shown in the following manner: since V ■ u — 0, 
the velocity field can be expressed as u = V x ip, where tp is the 
velocity vector potential with the gauge condition V • ip = 0. We 
then have u k {q) = iq p e kp iipi(q). Substituting this expression in 
Eq. (20), and using the fact that the velocity field is divergence- 
less, we get 



£i(k) = kj I B °_ 2 q Uj(k - q)i>i(q) dq 



J 



m 

Bp q 

m 2 



qi ui(k - q)ipi(q) dq. 



(21) 



For homogeneous and isotropic turbulence, the Uj{k — q)"4>i(q) 
correlation is proportional to S i (k). Since the first term in Eq. (21) 



is oc hi, it does not contribute to £ . Therefore, 



£i(k) 



Bom <ln 

m 2 



qi ui(k - q)ipi(q) dq. 



(22) 



Again, for homogeneous and isotropic turbulence, the 
Ui(k — q)tpi(q) correlation is oc 8 (k)g(\q\). One can 
then carry out the angular integral in Eq. (22) using 
J(q m qi/q 2 )(dQ/4n) = (l/3)<5 mi to get £i(x) = aB oi , 
where a is given by 



3?7 



(23) 



which is identical to the expressions obtained by Krause & Radler 
(1980); see also Radler & Rheinhardt (2007). (Note that when 
the Lorentz force becomes important the assumption of isotropy 
in the above derivation breaks down and the a-effect becomes 
anisotropic, as calculated below and detailed in Section 3.5.) 

Since we have already solved for the velocity field explicitly, 
we can now derive an expression for the mean emf in this approach. 
Substituting the velocity in terms of the forcing function, the mean 
emf in coordinate space is given by 



£i(x) 



f Bo g 

J m 2 



F jk (q) 



\vq 2 + (B„ ■ q) 2 lm 2 \ 



dq. 



(24) 



Here F jk the spectrum tensor for the force-field given by Eq. (15). 
Note that only the antisymmetric part of F jk contributes to £ , due 
to the presence of tij k on the RHS of the above equation. We can 
also write £ as 



£i(x) 



>o • q 



Fjk{q) 



(m 2 )(vq 2 ) 2 [i + Af 



dq, 



(25) 



where N = (Bo • q) 2 /(j)vq*) determines the importance of the 
Lorentz forces on the mean emf. It is to be noted that the limit of 
small Lorentz forces corresponds to taking N 1 above. 



3.1 Method A: express b in terms of u and then solve for £ 

In this approach we use the induction equation to solve for b in 
terms of u. Using Eq. (10) to express b in terms of u in Eq. (19) 
gives 



/" Bp- 

£lJfe J m 2 



£i(k) 
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- u 3 (k~ q)u k {q) dq. 



(20) 



3.2 Method B: compute £ directly 



In this approach, we directly compute £ = u x b by substituting 
u and b in terms of /, using Eqs (12) and (13). We then get 



£i(k) 



itijk J 



Bo q fj(k-q)f k (q) 



m 2 l(k - q)l(q) 



dq, 



(26) 
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where we have defined 7(g) = vq 2 + (Bo ■ q) 2 /r/q 2 . Substituting 
for the force correlation, the mean emf in coordinate space is given 
by 



B q 



F jk (q) 



dq, 



m 2 [vq 2 + {Bq ■ q ) 2 lm 2 f 
which is identical to Eq. (24) for £ obtained by Method A. 



(27) 



3.3 Method C: express u in terms of b and then solve for £ 

Note that one could also start from the momentum equation to com- 
pute £ . In this approach, we first solve for u in terms of b and the 
forcing function /, and then substitute for b in terms of / using 
Eq. (13). The difference from the earlier treatments will be an ad- 
ditional term containing an / x 6-like correlation, which turns out 
to be essential for calculating the £ correctly. Using Eq. (11) one 
can write 



Ui{k) = — [(ifc ■ Bo) bi(k) + fi(k)] . 

From Eq. (19) the mean emf can then be written as 

£i(k) 



(28) 



:fj(k - q)b k (q) dq 



1} 

(k-q)f 



v(k 



bj(k - q)b k (q) dq. 



(29) 



Here the first term involves the / x 6-like correlation. To elucidate 
the meaning of the second term it is useful to define the magnetic 
field as b = V x a, where a is the small scale magnetic vector 
potential in the Coulomb gauge (V • a — 0). Then, for isotropic 
small scale fields, following the approach in Method A, the second 
term in Eq. (29) gives a contribution to £ of the form qm-Bo, where 



qm 



1 r 

— ab. 

3v 



(30) 



So this contribution is proportional to the magnetic helicity of the 
small scale magnetic field (analogous to the helicity of the vector 
potential of the velocity field). 

If we substitute b in terms of / from Eq. (13) and then inte- 
grate over the delta function, the mean emf in coordinate space can 
be expressed as 



£i(x) = ieijfc 



Bo q 



Fjk(q) 



f Bp q 

J m 2 



(vq 2 ) 2 [1+N] 
Fjk(q) 



dq 



(31) 



N 



(vq 2 ) 2 [1 + N] 1 + N 



dq. 



The two terms on the RHS of the above equation have an interest- 
ing interpretation. As mentioned above, the limit of small Lorentz 
forces corresponds to taking N = (Bo ■ q) 2 /r/vq 4 <g; 1. In this 
limit the second integral vanishes while the first one [i.e. the f xb- 
like correlation in Eq. (31), which is really a (V -2 /) x b correla- 
tion] goes over to a kinematic a-effect. [One can see by comparing 
Eq. (25) and the first term in Eq. (31) that the two are identical in 
the limit N <g; 1]. In fact, this part of the a-effect can be obtained 
from Eqs (12) and (13) by neglecting the Lorentz force in the ex- 
pression for Wj(fe). In the following we refer to the contribution of 
the field-aligned component of this term divided by Bo, as the qf 
term, because it comes from the / x b-like correlation. 

As N is increased the contribution from the first term de- 
creases. In the same limit, the second term, which depends on the 
magnetic helicity, gains in importance. Since it has the opposite 
sign, it partially cancels the first term and further suppresses the 



total a-effect. This is reminiscent of the suppression of the kinetic 
alpha due to the addition of a magnetic alpha (proportional to heli- 
cal part of b) found in several closure models (Pouquet et al. 1976; 
Kleeorin & Ruzmaikin 1982; Gruzinov & Diamond 1994; Black- 
man & Field 2002; Brandenburg & Subramanian 2005b). 

When adding the two terms in Eq. (31), the mean emf turns 
out to be 



£i(x) 



[ Bp q 
J W 



FAl) 



[vq 2 + (Bo ■ q) 2 /m 2 ] 



dq, 



(32) 



which is identical to the expressions obtained in Methods A and B; 
see Eqs (24) and (27), respectively. 



3.4 Method D: compute £ from the d£ /dt — relation 

In recent years the so-called r-approximation has received in- 
creased attention (Kleeorin et al. 1996; Blackman & Field 2002; 
Radler et al. 2003; Brandenburg & Subramanian 2005a; Radler & 
Rheinhardt 2007). This involves invoking a closure whereby triple 
correlations which arise during the evaluation of d£/dt, are as- 
sumed to provide a damping term proportional to £ itself. In the 
present context there is no need to invoke a closure for the triple 
correlations, because these terms are small for low fluid and mag- 
netic Reynolds numbers. It turns out that the correct expression for 
£ can still be derived in the same framework, where one evaluates 
the d£ /dt expression. 

The expression for d£/dt is governed by two terms, u x b 

and u x b, where dots denote partial time differentiation. Of 
course, both ii an b vanish in the present case, but this is the result 
of a cancellation of driving and dissipating terms. In the present 
analysis both terms will be retained, because the dissipating term, 
which is related to the desired £ , can then just be written as the 
negative of the driving term. 

We perform the analysis in Fourier space and begin by defin- 
ing E(k, q) = u(k — q) x b(q). Note the required 



/ 



E(k,q)e lkx dkdq. 



(33) 



To calculate the time derivative d£ / dt, one needs to evaluate E = 
£k + Em, where 



E K (k, q) = u(k - q) x b(q) 



E M (k, q) = u(k - q) x b(q). 



(34) 
(35) 



For ii and b we restore the time derivatives in Eqs (10) and (11), 
and obtain 



E K = iq ■ Bo [u(k - q) x u(q)] - r/q E, 

E M = i(fe - q) ■ Bq [b(k - q) x b(q)] - v(k - q) 2 E 



(36) 



+/(fc - q) X b(q). 



(37) 



Since all time derivatives are negligible, we can simplify the RHS 
of the above equations by using Eqs (12) and (13) to express u and 
b in terms of the forcing function. Adding the two Eqs (36) and 
(37) yields 

■3 fu \iQ' Bo tijk Fjk(q) 



[ V q 2 + v(k-q) 2 ]E t = 8 3 (k) 
X m 2 +l{k-q) 



r/q 2 7(<7h(fc - q 
((k-q)-Bo) 2 
ri(k - q) 2 



(38) 
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where the function 7 was defined just below Eq. (26). The expres- 
sion in the squared brackets on the right hand side of the expression 
above exactly reduces to r\q 2 + v(k — q) 2 and so, in the steady state 
limit, d/dt — 0, we have 



E i (k,q) = 5 3 (k) 



iq ■ Bp £jjk Fjkjq) 

m 2 7(<?h(fc - q) ' 



(39) 



Using this expression in Eq. (33), and integrating over fc, we again 
recover the form of £ identical to Methods A, B and C. 

Thus, in this simple example where one can apply FOSA to 
both the induction and momentum equations, one gets identical ex- 
pression for £ in terms of the correlation properties of the forcing 
function /, in all the four methods. 

3.5 The nonlinear a-effect 

We now compute the nonlinear a-effect explicitly from the expres- 
sion of £ as obtained in the four methods discussed above. As has 
been mentioned earlier, only the antisymmetric part of Fjk con- 
tributes to £ in Eq. (15), so Eq. (32) takes the simple form 



£i{x) 



f Bp q 

J m 2 \ 



dq, (40) 



87rg 4 [vq* + (Bo ■ q) 2 /VQ 2 ] 

Contraction between the two e's and solving for a = £ ■ B /B 2 
leads to 



/ 



x(q) 



(B -q) 2 



dq 



W V [l + (Bo • q) 2 /3 2 ] 47R T 



(41) 



where we have introduced (3 2 = B 2 /(r\vq 2 ) and hats denote unit 
vectors. The solution involves an angular integral with respect to 
the cosine of the polar angle, /i — Bo ■ q, 



HP) 



so that 



fi 2 dfi 



+ /3V) 2 



1 



tan -1 /? 



1 + / 



2rju 2 



(42) 



(43) 



Note that for small values of (3 we have F((3) w 2/3 
the limit of large values of B and (3 we have F(f3) - 
so the expression of a reduces to 



* — f° 

^_ fW / 

4B 3 V v J 



^dq (for Bo 



00). 



4/3 2 /5. In 
V(2/3 3 ), 



(44) 



So, in the asymptotic limit of large Bo, we have a — > Bq . This 
is a well-known result that goes back to the pioneering works of 
Moffatt (1972) and Rudiger (1974); see also Rudiger & Kitchatinov 
(1993). Note that £ x Bq = 0, because the corresponding angular 
integral would be over the product of a sine and cosine term which 
vanishes. 

To illustrate further the dependence of a on Bo we need to 
adopt some form for the spectrum xil)- We assume that the forcing 
is at a particular wavenumber, go, and choose x(l) = H{5(q — go) 
where H f is the helicity of the forcing. Then the integration of the 
delta function simply gives 



a_ _ 3 / Bp V 
a ~ 2 VBcJ 



tan" 1 (Bo /.Be: 



Bo/ B c 



1 + B 2 /B 2 



where we have defined 

a = -Ht/{3r)v 2 ql), 



**Jwqo- 



(45) 



(46) 
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Figure 1. Variation of a, & F , and -a M with Bo from the analytical the- 
ory. Note the mutual approach of &p and -5m (asymptotic slope —2) to 
produce a lower, quenched value (asymptotic slope —3). 



If we express a — Qp + &m, where ap is computed from the 
(V~ 2 /) x b term and q m from the (V~ 2 f>) x b term in Eq. (31), 
we have 



-=3(|o\ 

ao \ B CI ) 



oto 

Q 



tarr^Bo/Bc 



= 3 



Bo_ 

B CI 



Bo/B cr 

3tan _1 (B /Bc 
2" 



3/2 + B 2 /B c 2 



Bo/ B c 



B 2 /B 2 



(47) 



(48) 



The hats on the as indicate that a special choice has been made to 
divide a up into different contributions. A different choice with- 
out hats that had been derived under the r approximation, will be 
discussed in Section 4.5. 

We plot in Fig. 1 the variation of a with Bo. This shows that 
a ~ ao for Bo < B cr and in the asymptotic limit a decreases 
oc Bq" 3 . This figure also shows the variations of of and olm with 
Bo as predicted from Eqs (47) and (48). Both decrease asymptoti- 
cally like Bq 2 because here, unlike in Eq. (45), the term in squared 
brackets remains constant. Their sum decreases as B,^ 3 . One could 
also define a kinetic &k from the (V~ 2 w) x u term in Eq. (25). In 
this case, for steady forcing we have a = &k = or + &m- 



3.6 Comparison with simulations 

Simulations allow us to alleviate some of the restrictions imposed 
by the analytical approach such as the limit of low fluid and mag- 
netic Reynolds numbers, but they also introduce additional restric- 
tions related for example to the degree of anisotropy. We adopt here 
a simple and frequently used steady and monochromatic forcing 
function that is related to an ABC flow, i.e. 

/ sin kfZ + cos k{ y \ 
f(x) — — L sin ki x + cos fa z ] , (49) 
* \ sin kty + cos kfX J 

where Af denotes the amplitude and k{ the wavenumber of the forc- 
ing function. This forcing function is isotropic with respect to the 
three coordinate directions, but not with respect to other directions. 
The helicity of this forcing function is H{ = k { A 2 . We use the 
PENCIL CODE, 1 which is a high-order finite-difference code (sixth 



http : / /www . nordita . dk/ software/pencil- code 
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0.001 0.010 0.100 1.000 10.000 

Bo / B C r 



Figure 2. Variation of a, d F , and — Am with Bo for the ABC flow at low 
fluid and magnetic Reynolds numbers (Re = R m = 10~ 4 ), compared 
with the analytic theory predicted for a fully isotropic flow. Note that the 
numerically determined values of a-p are smaller and those of —&m larger 
than the corresponding analytic values. However they still add up to satisfy 
the relation ctp + Am = a, predicted by analytic theory. In all cases the nu- 
merically determined values of &K agree with the numerically determined 
values of a. 



8 











1.00 








0.10 






Sf+Sm 


0.01 






0(K/(l+ofi m ). 
a=0.7 



0.1 1.0 10.0 100.0 1000.0 



Figure 3. Dependence of the a-effect on Rm for fixed field strength, Bo = 
u re f , for ABC-flow forcing at Re = 1. Note the agreement between a and 
op + &m as well as »k/(1 + aRm) for i? m < 30. 

mentum equation and hence the viscous time scale. However, since 
Re is small, no correction is necessary here. Again, a numerical 
resolution of 32 3 meshpoints was used except for Rm > 10, where 
we used 64 3 meshpoints. 



order in space and third order in time) for solving the compressible 
hydromagnetic equations. We adopt a box size of (2n) and take 
Ac — 10~ 4 , fef = 1, and determine a = £ z /Bo z as well as &k, 
&f, and &m, which are given respectively by the three integrals in 
Eqs (20) and (29). The result is shown in Fig. 2. For these runs a 
resolution of just 32 3 meshpoints is sufficient, as demonstrated by 
comparing with runs with 64 3 meshpoints. 

For Bo /Bct < 1 the resulting values of a agree in all cases 
perfectly with both &k and &f + &m- However, for Bo/B cr > 1 
the numerically determined values of a are smaller than those ex- 
pected theoretically using q = fef . As in the analytical theory, the 
quenching is explained by an uprise of — om. Note, however, that 
in the simulations this quantity attains a maximum at somewhat 
weaker field strength than in the analytical theory; cf. the dashed 
line and crosses in Fig. 2. We believe that this discrepancy is ex- 
plained by an insufficient degree of isotropy of the forcing function. 

Finally, it is interesting to address the question of the Reynolds 
number dependence of the quenching behaviour. In Fig. 3 we show 
the results for a, &k, &f, and —&m- Since the velocity is of the 
order of u re f = Af/(vkf ), we have defined the fluid and magnetic 
Reynolds numbers as Re = u le {/{vki) and R m = u Te i/(r}kt), 
respectively. For all runs we have assumed Re = 1 and Bo = u re f . 

Again, we see quite clearly the approach of —&m toward dp, 
so as to make their sum diminish toward a with increasing values 
of Rm - For Rm < 1 the numerical data agree well with the analytic 
ones, whilst for R m > 1 the numerical values for all alphas lie be- 
low the analytic ones (not shown here). In particular, for Rm > 1 
the value of qk, based on the integral in Eq. (20), begins to exceed 
the value of a. This apparently signifies the break-down of FOSA. 
However, one may expect that the relevant inverse time scales or 
rates are no longer governed by just the resistive rate, ~ rjkj, but 
also by a dynamical rate, ~ w re f fef . This leads to a correction factor, 
1/(1 + aRm), where a w 1 is an empirically determined coeffi- 
cient quantifying the importance of this correction. In Fig. 3 we 
show that both qf + &m as well as &k/{1 + aRm) with a = 0.7 
are close to a for R m < 30. Note that no correction is necessary 
for Of or &m, because these quantities are determined by the mo- 



4 TIME-DEPENDENT FORCING 

We now consider the case when / depends on time, but is neverthe- 
less statistically stationary. In that case, both b and ii are finite, and 

hence also u x b and lixi can in general be finite, even though 
their sum might vanish in the statistically steady state. Later we spe- 
cialize to one case of particular interest, when the correlation time 
of the forcing is small. This was the case, for example, in the sim- 
ulations of Brandenburg (2001) and Brandenburg & Subramanian 
(2005b). As in Section 3, we will assume here small magnetic and 
fluid Reynolds numbers and neglect the nonlinear terms in the in- 
duction and momentum equations, but retain the time dependence. 
We will also now take a Fourier transform in time and define 



u(k, co) 



u(x,t)e- lk - x+ilot dxdt 



(277)4 

which satisfies the inverse relation 
u(x,t) = / u(k,Lj)e k ' x ~ 1 ^ t dkdu>. 



In Fourier space, Eqs (10) and (11) become 

(— iw + r)k 2 )bi(k,u>) = (ife • Bo) Ui{k, to), 

{— iw + vk 2 )ui{k, u>) = (ife ■ Bo) b»(fe, u>) + f%(k, u>). 



(50) 



(51) 



(52) 
(53) 



In order to simplify the writing of the equations below, it is conve- 
nient to define complex frequencies 



-iuj + r]k , r„ 



-iu) + vk 



(54) 



As before, we can solve equations (52) and (53) simultaneously to 
express u and b completely in terms of /, 



Hi(k, lo) — 
h(k,uj) = 



fi{k,L 



r„ + {B () -k) 2 /r v ' 

fi(k, lo) ifc ■ Bp 

r„ + (Bo ■ k)yr v r v ' 



(55) 
(56) 
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We can use these solutions to calculate £ . For getting an explicit 
expression, we also need the force correlation function in Fourier 
space. For isotropic, homogeneous and statistically stationary forc- 
ing, this is given by 



fj(p,v) fk(q,oj') =5 (p + q)S(u + J)F jk (q,u), 



(57) 



where we can still take F jk to be of the form given by Eq. (15), 
with spectral functions now changed to say, a frequency dependent 
&(k,co) mdx(k,u;). 

Note that in the limit where the correlation time of the forcing 
function is short (delta-correlated in time) the Fourier space spec- 
tral function Fjk is nearly independent of the frequency w. How- 
ever if one evaluates the helicity of the forcing, one gets 



X(k,u) dkduj = f-(Vxf) = H t , 



(58) 



where the k integral is from to oo, while the uo integration goes 
from —oo to +00. For \ independent of cu this would be infinite. So 
we still keep a spectral dependence and write x(q, uj) = x(<7)g(w), 
where g is an even function of u>, satisfying J g(uj)du> = 1. [The 
property that g is even is a consequence of the forcing function 
being real; see Eq. (7.44) of (Moffatt 1978).] For a forcing with say 
a correlation time t, g(u) will be nearly constant for ujt ~ 1 and 
decay at large u. In the limit of small r the range for which g(ui) 
is nearly constant will be very large. We will need only g(0) ~ r 
in most of what follows. Note that the other extreme limit of steady 
forcing corresponds to taking g(cu) — > 5(u>). The mean emf can be 
written as 



£i(x, t) = £ij k Ujb k 



= j £i(k, 



*dfc dfi, 



where the Fourier transform £ is given by 



£i(k, = Ujk \ Uj(k — q,fl — Lu)b k (q,u) dqdco. 



(59) 



(60) 



We now turn to the calculation of the nonlinear mean emf and 
the resulting nonlinear a-effect. We focus on Method A, the di- 
rect Method B and also Method C to illustrate the similarities and 
differences from the case when the time evolution is neglected. We 
also discuss in detail the result of applying a r-approximation type 
method in the subsequent section. 



4.1 Computing £ from the induction equation 

As before, we first start from the induction equation to solve for b 
in terms of u using Eq. (52), so 



c 11 \ ■ [ T3 Uj(k-q,Q.~u))u k (q,u>) 

£i(k) = ie i3fc / B q — : ■ dqdco. (61) 

We can then express u in terms of / using Eq. (55). Substituting 
from Eq. (57) for the force correlation in a time dependent flow, the 
mean emf in coordinate space is then given by 



B q 



Fjk{q,u) 



\Tu + (Bo ■ q) 2 /T v 



-r dq dto . 



(62) 



As usual, we define a = £ ■ Bq/Bq. Since only the antisymmetric 
part of Fjk contributes in the above integral, we have, 



q = _ [ (Bo-q) 2 
J 4nq' 



? 4 1 r„r„ + (Bo • q) 2 1 2 
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In the following we refer to this expression for a also as a K , since 
it is seen to arise purely from the (V _2 ux w)-type velocity correla- 
tion, generalized to the time dependent case. Note that the denom- 
inator of Eq. (63) is even in cu and so is the spectral function g(u>). 
Therefore the term in the above integral, which has an iui in the 
numerator, vanishes on integration over w (by symmetry). So the 
mean emf is a real quantity as it should be. Before evaluating the 
above integral explicitly, let us ask if we get the same expression 
for Eq. (62) using Methods B and C, even in the time-dependent 
case. 



4.2 Computing £ directly 



Let us directly compute £ = u x b by substituting u and b 
in terms of /, using Eqs (55) and (56). We also substitute from 
Eq. (57) for the force correlation in a time dependent flow. The 
mean emf in coordinate space is then given by Eq. (62), so we do 
not repeat it here. 

4.3 £ from the momentum equation 

As before we start from the momentum equation, solve for u in 
terms of b and the forcing function /, and then substitute for b 
in terms of / using Eq. (56). We particularly wish to examine if 
the (V~ 2 /) x £>-like correlation is essential for calculating the £ 
correctly, even for the time dependent case. Using Eq. (53) one can 
write 



iii(fc,w) = [(ik-Bo)bi(k,u) + fi(k,u)] . 
From Eq. (60) the mean emf can then be written as 



(64) 



£i(k,n) 



eijkl Mk / - q ' n : uj ]t (q :Z ) dqdc 



-i(fi-aj)+^(fc-q) 2 

f Bo • (fc - g) 

+KiSk J - lj) + u{k - q) 2 

x bj(k — q, Q, — to)bk(q, u>) dq dco. (65) 

Here the first term involves the (V~ 2 /) x b-like correlation, gen- 
eralized to the time-dependent case. Substituting b in terms of / 
from Eq. (56) and integrating over the delta functions in wavenum- 
bers and frequencies, the mean emf in coordinate space can be ex- 
pressed as 



£i(x) 



. [Bo- 



q F jk (q,oj) 



1^27 k 



r„| 2 [i + n] 

Bp-q F jk (q,io) 



dq duo 

N* 



(66) 



r„ irj 2 \i + N\ i + n* 



dqdco. 



Here we have defined N = (Bo ■ q) 2 /T n Y v . Now the limit of 
small Lorentz forces corresponds to taking \N\ <t; 1. Again in 
this limit the second integral vanishes while the first one [i.e. the 
generalized (V~ 2 /) x 6-like correlation in Eq. (66), goes over 
to a kinematic a-effect. In fact, this part of the a-effect can be 
obtained from Eqs (55) and (56) by neglecting the Lorentz force 
in the expression for Ui(k, oj). 

On adding the two terms in Eq. (66), the mean emf turns out to 
be identical to Eq. (62), obtained when starting from the induction 
equation or in the direct Method. Therefore, one could either com- 
pute the mean emf starting from the induction equation, or directly, 
or from the momentum equation, by the addition of a generalized 
(V~ 2 /) x 6-like correlation and a purely magnetic correlation. 
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We can again define ap and &m for the time dependent forc- 
ing, from the first and the second terms on the r.h.s. of Eq. (66) 
respectively. We have 



QF 



QM = 



/ 



(Bp ■ g) 2 x (q) gH [r;r; + (go ■ q ) 2 } dqdm 
4tt ? 4 r; | r„r„ + (b ■ q) 2 \ 2 



,(67) 



gg (go ■ q) A X {q) gH dqd^ 

47r g 4 r; ] r„r„ + (B ■ q) 2 I 2 



(68) 



In terms of the total a, we have a — a F + a M - It is now explicitly 
apparent that in the time dependent case, a = &k — &f + Am in 
agreement with what is obtained for steady forcing. Further, in the 
limit of steady forcing where, g(oj) — > S(oj), the above general- 
ized expressions reduce to the a's obtained from Eqs (25) and (31) 
respectively. 



4.4 The nonlinear a-effect 

Let us return to the explicit computation of the nonlinear a-effect 
for the delta-correlated flow. Solving for a — £ ■ Bo/B 2 from 
Eq. (63) leads to 



n= - I X ( q )(B -q) 2 ( V q 2 )I(q) d " 



4nq 2 



(69) 



|r„r„ + (B () ■ q ) 2 



(70) 



where / is the integral over oj given by 
g(oj)doj 

|r„r„ + (Bo ■ q) 2 \ 2 = 

where the latter expression for I obtains in the limit of a delta- 
correlated forcing, where g(oj) = r is almost constant through the 
range where the rest of the integrand contributes significantly. We 
now focus on the special case of v — r\ for the explicit evaluation 
of the above integral. Note that calculating / is then straightfor- 
ward but tedious. We briefly outline the steps and then quote the 
result. First the denominator of Eq. (70) can be expanded and then 
factorized to give 



|rvr„ + (B • qf\ 2 = (oj + z)(oj + z*)(oj-z)(oj-z*), 



(71) 



where z — (Bo ■ q) + \vq 2 . Then the integral can be rewritten as 



rdw 



2(2* 2 - Z 2 )\z\ 



oj — Z* OJ + z* 



- + - 



Z OJ + z 



(72) 



Grouping the term having oj — z* with the one having oj — z and 
likewise the two terms with uj + z* and cu + z, we get 



F 

J — c 



T&LO 



2(z*+z)\z\ 2 
-oj + 2Bo ■ q 



(73) 



+ ■ 



oj + 2B ■ q 



(oj - Bo ■ q) 2 + v 2 q 4 (oj + Bo ■ q) 2 + v 2 q 4 

where the last expression is explicitly real. A change of variables 
allows the above integral to be done easily to give 2 



/ = 



(74) 



2vq 2 [(Bo-q) 2 + v 2 q A Y 

Substituting / into Eq. (69), and carrying out the angular integral 
over fi — Bo ■ q then gives 



2 We thank K.-H. Radler for pointing out to us that this result can be gen- 
eralized for r] ^ u, to give I = ttt/{(t] + v)q 2 [(Bo ■ q) 2 + rjuq 4 ]}. 



poo 

2 Jo » 



X ^G(P)dq, 



where f3 = Bo jvq for the v = r/ case, and 



G(P) = ^ 



1 - 



tan" 1 /; 



(75) 



(76) 



So for the time dependent delta-correlated flow, in the asymptotic 
limit of large Bo, we have a — > B^ 2 . If we further assume that 
the forcing is at a particular wavenumber, q , and choose x(<?) = 
HfS(q — go), we have 



a 
eta 



(.Bcr) 



tan-^Bo/Bcr 



Bo /B c 



where we have now defined 

7T TH{ 



ao 



6 q 2 B 



2 • 



vqo- 



(77) 



(78) 



4.5 Comparison with r-approximation 

Note that in the case of time dependent forcing, considered for 
example in the simulations of Brandenburg (2001) and Branden- 
burg & Subramanian (2005b), both b and u are finite, and hence 

also u x b and it x b can in general be finite, even though their 
sum would vanish in the statistically steady state. So taking a time 
derivative of £ and then examining the stationary situation, could 
break the degeneracy between a K and q f in the kinematic limit, 
and also lead to novel insights. Indeed, if one were not able to solve 
for u and b explicitly in terms of the forcing this would be the prac- 
tical route to follow. 

We now examine the time-dependent case in a manner analo- 
gous to the so-called r-approximation. As mentioned earlier, the r- 
approximation closures involve invoking a closure whereby triple 
correlations which arise during the evaluation of d£/dt, are as- 
sumed to provide a damping term proportional to £ itself. In the 
present context there is no need to invoke a closure for the triple 
correlations, because these terms are small for low fluid and mag- 
netic Reynolds numbers. It turns out that the correct expression for 
£ can still be derived in the same framework, where one evalu- 
ates the d£ /dt expression. Further, it allows us to define a new set 
of a's for the time dependent forcing, an, cum and ckf, i.e. with- 
out hat, which respectively incorporate the kinetic (u, u), magnetic 
(b, b) and the force-field (f,b) correlations (see below). These 
a's have properties very similar to those which arise in the r- 
approximation closure for large R m systems. 

For example, we showed above that for steady forcing, the 
V -2 / x b correlation is non-zero and essential to calculate the £ 
correctly in Method C. We also showed above that such a correla- 
tion is important to include even for a time-dependent forcing, if 
one starts with the explicit solution of the momentum equation for 
u. On the other hand, in a r-approximation type approach one first 
evaluates the d£ /dt expression, which involves a ii x b type cor- 
relation, instead of solving first for u. It is then interesting to ask, 
is the corresponding (/, b) correlation, or the ap term defined be- 
low, which arises in the evaluation of u x b, still non-zero for the 
time-dependent, delta-correlated forcing? Or does it vanish in the 
r-approximation type approach, as assumed in earlier work (Black- 
man & Field 2002; Radler et al. 2003; Brandenburg & Subramanian 
2005a)? In particular, does one then recover a = qk + «m, a re- 
lation which one obtains in the r-approximation at large R m 7 We 
examine these issues in detail below. 
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We write the time derivative of the emf as £ = £ k + £n 



where £k = u x b and £m = u x £>. From the induction equa- 
tion for b and the momentum equation for u, we now have 



£ K = wxfo = uxB - Vu + u x 77V 2 6, 



(79) 



£ M = ux6 = /xb + BoV6x6 + iA7 2 u x 6, (80) 

where the perturbed pressure term vanishes for divergence free 
forcing. We can evaluate each of these terms in Fourier space, since 
we have the Fourier space solutions for both u and b completely in 
terms of /. 

In the time dependent case, we define by analogy to Eq. (33), 



E(k, q, Q, uj) = u(k — q,Q — u>) x b(q,u)), (81) 
so that the emf in coordinate space is, 



£ = J E{k,q,Q,w) e ikx - iClt dkdqd0.du. 
We also define, analogous to Eqs (34) and (35), 



(82) 



E^\k,q,n,u) = u(k-q,n-Lu)x [-iw6(g,w)], (83) 
E$(k,q,Q,w) = [-i(fi-w)u(fc-qr,fi-w)] x 6(g,w).(84) 

Note that + E$ = —iQE. Using the induction and momen- 
tum equations, we have explicitly, 



e£> = iq ■ B [u(k - q,Q - w) x u(q,u))] - rrf E (85) 



E 



(t) 



i(k — q) ■ Bo [b(k — q,Q — w) x b(q, w) 



Substituting -it and 6 in terms of / from Eq. (55) and (56), and 
integrating over the delta functions in wavenumbers and frequen- 
cies which arises in taking the (/, /) correlations, we then have in 
the coordinate space, 



qk 



/ 



QM 



{Bp-qf X (q) ffHlFj 2 dqdcj 
4nq 4 | r„r, + (Bo • q) 2 \ 2 ( V + v)q 2 ' 

BtjBo ■ q) 4 X (q) g(u) dgdco 
4<Kq 4 \ T„T n + (B ■ q) 2 \ 2 (n + v)q 2 ' 

(Bo ■ q) 2 X (q) <?H [rjr; + (b ■ g ) 2 ] d g du 



4nq 4 j r„r^ + (Bo ■ q) 2 \ 2 (r) + v)q 2 
On adding Eqs (94), (95) & (96), the expression for a is 

\2 



[(Bp ■ q) 2 x(g ) g( U ) I 
J 4nq 4 (r) + v)q 2 



where 



I a = 



(B ■ g) 2 + r»r; + (go ■ g) 2 



+ f(k-q,n-iu)xb(q,uj)-v(k-qyE (86) 
We add Eq. (85) and Eq. (86), to write 

E 

—iQE H = iq- B ®(u,u) +i(k - q) ■ B $(b,b) 

T~cff 

+*(/,&), (87) 
where we have defined, for any pair of vector fields f 1 and / 2 , 

*(/i,/ 2 ) = f\(k-q,n-u) X f 2 (q,w), (88) 

and = r\q 2 + v(k — q) 2 . We note in passing that the above 
equation is similar to the corresponding equation which obtains un- 
der r-approximation in the large Reynolds number case, except that 
T ef f would then correspond to a relaxation time for triple correla- 
tions (cf. Brandenburg & Subramanian 2005a). So we have, 

E = r* ff [iq ■ B $(u, u) + i(fe - q) ■ B $(6, b) + $(/, b)] .(89) 

where we define t* s = r e ff/[l — iSl T e s]. 

Letusdefinea = (£-Bq)/B 2 as before. Then in coordinate space, 



(94) 
(95) 
(96) 

(97) 
(98) 

b 2 J I r„r„ + (Bo ■ q) 2 | 2 ~" (99) 

It is apparent from Eq. (99) that the expression for a turns out to 
be the same as in Eq. (63). Therefore, the r-approximation type 
treatment also gives the same a as the other methods. 

It is interesting to consider what happens to the various a's 
defined in Eqs (94)-(97), in the limit of steady forcing, where 
g(cu) — > 8(lj). It is straightforward to check that in this steady 
forcing limit a defined in Eq. (99) goes over exactly to the total 
a = a' s ' given by the steady state expression in Eq. (41) of Sec- 
tion 3.5. Also, in the steady forcing limit, we get 

Qk -> lv/(7l + v)}aK = {v/(r 1 + v)W S) , (100) 
«m — » [v/{r) + v)]&m, (101) 
a F -» [v/(r) + v)]&F. (102) 

In this limit one has therefore 

Qm + qf = [v/(;n + f)](a M + or) = W/(v + v)]a {s) (103) 
and so, once again, 



r„r„ + (B ■ q) 2 | 2 

The numerator of this integrand can be simplified to give 
T* [r/q 2 + vq 2 \ such that, 

(Bo • q f x (q) g(u>) r; d(?d ^ 



Qk + Qm + a-p = a 



(S) 



(104) 



Q = Qk + Qm + Qf, 
where 



(90) 



a K = J iq ■ Bo $(«,«) • B T* g e ik x - int dkdqdndoj, (91) 
a M = iq- Bo $(6, b) ■ B r c * ff e ifeiK " int dfc dq dfi dco, (92) 



as expected. It should be emphasized, however, that for a general 
time dependent forcing, there is no simple relation of the form 
given by the expressions (100)-(102). 

Let us consider the case of delta-correlated forcing now in 
more detail. It is of interest to check if the qf term contributes 
in the r-approximation type closures, as it does in the time- 
independent case. We have from Eq. (96), 

' x(q) (Bo-4) 2 



Of = 



■Hq) 



dq 



(105) 



a F = J iq- B $(/, b) ■ B r* ff e'^^'Mfc dq dQ dco (93) 

correspond respectively to the terms containing the (u, it), (b, b) 
and (/, b) correlations on the RHS of Eq. (89). 



(r) + v)q 2 ryH/ 4nq 2 ' 
where If is the integral over uj given by 

If = I ,„ v „ v r=r- — -)z i „ q(cu)duj. 

J ITvTr, + (Bo ■ q) 2 | 2 yK ' 

Here we can simplify T*r* + (B • q) 2 — —u> 2 + vr\q 4 + (Bo ■ 
q) 2 + iu(r)q 2 + vq 2 ). The integral over the term odd in w again 
vanishes, leaving again a real 7f, 



(106) 
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If 



I 



-lj 2 + urjq 4 + (Bp ■ q) 2 
\T„T V + {B ■ q)2|2 



ff(w)dw. 



(107) 



Let us focus on the case r\ = v as before. We rewrite the numerator 
using the identity —uj 2 + vriq 4 + (Bo -q) 2 = — (lj + z)(u) — z*) + 
uj(z — z*), so we have 

* — (lj + z)(lo — z*) + lj(z - 



If 



g(uj) duj. 



(108) 



r„r v + (B - q ) 2 \ 2 

The second term in the numerator of Eq. (108) does not contribute 
to the integral, since it is odd in lj, while the denominator is even. 
To simplify the integral further we use the identity in Eq. (71) for 
its denominator, giving 

g(uj)d(j 



If 



F 

J — C 



f 



(u) + z*)(uj — z) 

g(uj) <lu f 1 

z + z* 

g(uj) du (u 

z + z* (u — x) 



(109) 



: U + Z* 

- x) + iy 



(lj + x) + iy 



+ y 2 (u + x) 2 +y 2 



where we have defined x = Bo ■ q and y = vq , which are the 
real and imaginary parts respectively of z. Now changing variables 
to u — lj — x in the first term and u — lj + x in the second term 
we see that 



u + iy g(u + x) — g(u — x) 



z + z* 



u 2 + y 2 



du = 0. 



(110) 



Note that If — > in the limit of a delta-correlated forcing, where 
g(uj) = r is almost constant through the range where the rest of the 
integrand contributes significantly; that is g(u+x) — g(u—x) = r 
where the integrand contributes significantly, while g(u + x) — » 
and g(u — x) — » at large u. This can also be checked by doing 
the integral for 7f numerically. So, interestingly, o?f = 0. Thus, for 
a forcing which is random and delta-correlated in time, there is no 
contribution from the / x b type correlation! Thus, a is the sum of 
just two terms, a kinetic and a magnetic contribution which can be 
shown explicitly as follows. 

We note that Eq. (94) can be expressed as 



qk 
where 
Ik = 



4tt(77 + v)Bl 



(Bo ■ qf # I K (q) dq, 



LJ 2 + V 2 q 4 



rYr„ + (Bo ■ q) 



duj 



f LJ 2 + V 2 q 4 + (Bo-q 

j I r„r„ + (Bo 



) 2 - (So • q) 



duj. (112) 



q) 2 \ 2 

As before, we focus on the case r\ — v when the numerator can be 
simplified as lj 2 + ifq 4 + (B ■ q) 2 - (B ■ q) 2 = (lj + z) (lj - 
z*) - lj(z + z*) - (Bo ■ q) 2 . So we have 

duj. (113) 



Ik = r 



(lj + z)(l 



') - lj(z + z") - (Bo ■ qY 



T V V V + (Bo ■ q) 2 



It is to be noted that the second term in the squared bracket in 
Eq. (113) does not contribute to the integral. Using the identity in 
Eq. (71) for its denominator, we have, 



Ik = 



doj 



(lj — z)(lj — z*) 



F J(lj + z)(lj 



du) 



+ z*)(lj - z)(lj - z*) 



TYT 



(Bo ■ qf 



vq 2 2vq 2 (B ■ q) 2 + v 2 q A ' 



(114) 



Substituting Ik into Eq. (Ill) and carrying out the angular integral 
over fi — Bo ■ q then gives 



n 

a K = --r 



X(q) M + \rBl 



v 2 q 4 



where (3 = Bo/uq for the v = rj case, and 



1 _ J_ , i 

3 I3 2 ' 



tan _1 /3 



A similar analysis for Eq. (95) yields, 

S^oJ {B °- q)2] f lM{q) ^ 



QM 

where 

Im = 



(115) 



(116) 



(117) 



(Bo • q) 2 du 



V V 



V + (B - q) 
(Bo ■ qf 



2vq 2 (Bo- q) 2 + v 2 q 4 ' 
Carrying out the angular integral as earlier gives, 

x(q) 



71 d2 
ttM = -T-Do 



H(J3) dq. 



(118) 



(119) 



If we further assume that the forcing is at a particular wavenumber, 
go, and choose \(q) = H{8(q — go), we then have 



ao 

ttM 



1 3 

2 + 2 



(b c J 



1 




-2 " 




\B CI ) 





Bo / B cr 
tan" 1 (B /B c 

Bo /B cr 



(120) 



(121) 



where ao and B cr were defined in Eq. (78). It is explicitly apparent 
that a = (*k + qm, in agreement with Eq. (77). The result is plot- 
ted in Fig. 4. Note also that qf = 0, as was assumed in the minimal 
r-approximation (MTA) type calculations for large fluid and mag- 
netic Reynolds numbers (Brandenburg & Subramanian 2005a). 
It is interesting to note that in the limit Bo /B cr 3> 1, qk — » 
(111) +a /2 + O(B - 2 ) and a M -> -a /2 + 0(B^ 2 ) and so the total 
q = qk + cvm — > as Bq 2 . This is reminiscent of the kinetic and 
magnetic q's nearly cancelling to leave a small residual a-effect in 
EDQNM or MTA type closures. It is also interesting to consider 
the limit when Bo / B cr <C 1. In this limit qk — * Qo and olm — * 0, 
and so the net a-effect is just the kinetic contribution. Finally, for 
any B we note that (qk — aM)/ao = 1. 



4.6 Comparison with simulations 

It is appropriate to compare with the simulations of Brandenburg & 
Subramanian (2005b). We have produced additional results for low 
fluid and magnetic Reynolds numbers (Re = B m ~ 2 x 10~ 2 ). 
The forcing consists of helical waves with average wavenumber 
fcf/fci = 1.5. The resulting values of a are fluctuating strongly, 
so it is important to average them in time. Instead of calculat- 
ing the full integral expressions, we estimate the contributions to 
a from the formulae qk = —2t(u x u z , v ), oh = 2T{b x b ZiV ), 
q p = T if x b) ■ B /B 2 , where B = (0, B , 0) is the imposed 
field, and r _1 = (v + r/)k 2 . 

The result is shown in Fig. 5 and compared with the results 
of the previous section. In all cases the resulting values of of are 
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Figure 4. Variation of a, a K , and -a M with Bo from the analytical theory 
using a delta correlated forcing. Note that a = «k + «m and ap = 0. 
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Figure 5. Variation of a, c*k, and — qm with Bq for a random flow at low 
fluid and magnetic Reynolds numbers (Re = ij m ss 2 X 10 — 2 ), compared 
with the analytic theory predicted for a fully isotropic flow. Note that the 
numerically determined values of q?k are smaller and those of -om larger 
than the corresponding analytic values, which is similar to the results for 
the ABC-flow forcing. 



negligibly small and will not be considered further. Like in the case 
of the ABC flow, the numerically estimated values of a are smaller 
than the analytic ones. This might be explicable if for some reason 
the relevant normalization in terms of «o were to depend on Bo- 
Alternatively, the discrepancy might be due to us using only simpli- 
fied expressions instead of the full integral expressions. However, 
the important point is that the main contribution to the quenching 
comes from the growing contribution of -om such that qk + «m 
is quenched to values much smaller than «k. The corresponding 
results in the case of larger R m are given by Brandenburg & Sub- 
ramanian (2005b). 



5 DISCUSSION 

We have considered here the nonlinear a-effect in the limit of small 
magnetic and fluid Reynolds numbers, for both steady and time- 
dependent (both for general and delta-correlated) forcings. In the 
limit of low R m and Re, one can neglect terms nonlinear in the 
fluctuating fields and hence explicitly solve for the small scale mag- 



netic and velocity fields using double FOSA. We can then calculate 
the a-effect in several different ways. 

For both steady and time dependent forcings, one gets similar 
results, provided one starts from the explicit solutions to the induc- 
tion and momentum equations. Lets us begin with a summary of 
the results for the steady forcing case: To begin with we follow in 
Method A, the traditional route of solving the induction equation 
for b in terms of u, and then calculating the a-effect. For statisti- 
cally isotropic velocity fields, this gives a dependent on the helicity 
of the velocity potential, as already known from previous work. In 
addition since we have an explicit solution for u in terms of /, 
one can relate a directly to the helical part of the force correlation. 
In Method B we solved for u and b in terms of the forcing func- 
tion / and computed £ directly. This would correspond to what is 
done when the a-effect is determined from simulations. However, 
in general this cannot be done analytically unless one can solve 
for the small scale velocity and magnetic field explicitly. We get a 
identical to that obtained in Method A. 

More interesting is Method C, where one takes the momentum 
equation as the starting point, instead of the induction equation. In 
the limit of small fluid Reynolds numbers one can solve for u in 
terms of b and hence compute £ . This necessarily involves also the 
(V~ 2 /) x b correlation, between the forcing and the small-scale 
magnetic field, in addition to the (V _2 j) • b (or a ■ b) correlation 
arising from the Lorentz force. This second term depends on the he- 
licity of the small scale magnetic fields. When the Lorentz force is 
small, the first term contributes to a in a manner closely related to 
the usual kinematic alpha-effect, while the second term contributes 
negligibly. Interestingly, as the Lorentz force gains in importance 
the first term is suppressed, while the second term (which has an 
opposite sign) gains in importance and cancels the first term, to 
further suppress the total a-effect (Fig. 1). This is similar to the 
suppression of the kinetic alpha due to the addition of a magnetic 
alpha (proportional to the helical part of b) found in several clo- 
sure models (Pouquet et al. 1976; Kleeorin & Ruzmaikin 1982; 
Gruzinov & Diamond 1994; Blackman & Field 2002; Branden- 
burg & Subramanian 2005b). When one combines the two terms, 
the resulting a-effect is identical to that obtained from the induc- 
tion equation (Method A) using FOSA or the direct computation of 
Method C. However, it also highlights the fact that in this steady 
case the (V -2 /) x b-type correlation does not vanish, and that 
there is no tendency for this term to balance the viscous term, as 
one might have expected. 

Finally, the results of Method D show that the formalism used 
in the r-approximation lead to results that are equivalent to the 
usual approach taken in the first order smoothing approximation. 
However, this requires that the detailed spectral dependence of the 
diffusion operator be retained until the point where the steady state 
assumption is made. The resulting equations are solved for the 
spectral electromotive force of the form E(k, q) in Eq. (39). Other- 
wise, one would not recover the correct low conductivity limit, as 
shown by Radler & Rheinhardt (2007). We emphasize that through- 
out this paper we have understood the term r-approximation only 
in this more generalized sense. 

The above results are also obtained for statistically stationary 
but time-dependent forcing. Specifically, we showed that even in 
the time dependent case, one gets identical results for the mean emf 
if a is computed directly (as in Method B), or from the momentum 
equation, by the addition of a generalized (V~ 2 /) x fo-like corre- 
lation and a purely magnetic correlation. The explicit form of the 
a-effect differs between delta-correlated and steady forcing cases. 
In particular, in the limit of large Bo, and when v = r\, we have 
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Table 1. Summary of the results obtained from FOSA and r-approximation 
type analysis for steady and random forcings. 
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a oc B 2 for delta-correlated forcing, in contrast to a oc B 3 , for 
the case of a steady forcing. The former result has already been ob- 
tained by Field et al. (1999) and Rogachevskii & Kleeorin (2000), 
both of whom assumed the force-field correlation to vanish. How- 
ever, their result was derived under the assumption of large fluid 
and magnetic Reynolds numbers. 

The major difference between the time dependent and steady 
forcing cases arise when one follows Method D, the formalism used 
in the r-approximation type approaches to computing a in large 
Rm systems. We recall that in this approach one starts by evaluating 
the time derivative of the emf, and then look at the stationary limit. 
We showed that in this limit the a-effect can be naturally written 
as the sum of 3 terms, a = «k + cum + cuf, for a general time- 
dependent forcing. Here ok and ccm are the kinetic and magnetic 
contributions to a corresponding respectively to the terms contain- 
ing the (u, u) and (b, b) [see Eqs (91) and (92)], while qf incorpo- 
rates the (/, b) correlation; see Eq. (93). Interestingly, we showed 
that «f = in the approach of Method D, for delta-correlated 
forcing, and therefore a = qk + Qm, just the sum of a kinetic 
and magnetic terms. We also computed ckk and om explicitly for 
the case rj = v. In the kinematic limit, am — > 0, while B^ax. 
While in the limit B /B„ > 1, a K -> +a /2 + 0(Bq 2 ) and 
ttM — » — Qo/2 + 0(Bq 2 ), so that the total a — «k + ckm — > as 
Bq 2 . This is reminiscent of the kinetic and magnetic a's nearly 
cancelling to leave a small residual a-effect in EDQNM or r- 
approximation type closures. The results from employing FOSA 
and r-approximation type analysis is summarized in Table 1 both 
for steady and random forcings. 

As far as the low Reynolds number case is concerned, our an- 
alytic solutions demonstrate quite clearly that one can look at the 
nonlinear a-effect in several equivalent ways. On the one hand, 
one can express a completely in terms of the helical properties 
of the velocity field (Method A) as advocated by Proctor (2003) 
and Radler & Rheinhardt (2007). At the same time a can be natu- 
rally expressed as a sum of a suppressed kinetic part (first term in 
Method C) and an oppositely signed magnetic part proportional to 
the helical part of b (second term in Method C). 

Method D applied to the delta-correlated forcing is particu- 
larly revealing. As here one can explicitly write a — qk + «m, or 
as the sum of a kinetic qk, which dominates in the linear regime, 
and a magnetic cum, which gains in importance as the field be- 
comes stronger, and cancels qk to suppress the net a-effect. This 
is similar to the approach that arises from closure models like 
EDQNM (Pouquet et al. 1976) and the r-approximation (Klee- 
orin et al. 1996; Blackman & Field 2000; Brandenburg & Subra- 
manian 2005a) or the quasilinear models (Gruzinov & Diamond 
1994). In all these cases the nonlinear a-effect, for large R m , is 
the sum of a kinetic part and an oppositely signed magnetic part. 
As noted above, the kinetic a-effect is itself suppressed, as seen in 
Fig. 4, but this happens only for Bo > B cr , and the suppression is 
milder than the the strong suppression of the total a-effect. This is 



also borne out in the simulations of Brandenburg & Subramanian 
(2005b), where the kinetic part of the a-effect is suppressed in a 
manner that is independent of the magnetic Reynolds number, even 
though the total a is catastrophically suppressed. Finally, we also 
have shown that qf defined naturally in the approach of Method D, 
vanishes for delta-correlated forcing, as was assumed in derivations 
of the a-effect in r-approximation type closures. 

In the special case of periodic domains it is now clear that 
for forced turbulence at large B m the steady state a-effect is 
catastrophically quenched (Cattaneo & Hughes 1996; Brandenburg 
2001). However, the physical cause of this phenomenon was long 
controversial. Is it because Lorentz forces cause a suppression of 
Lagrangian chaos (Cattaneo et al. 1996) or is it due to a nonlinear 
addition to a due to helical parts of the small scale magnetic field, 
as is argued here? The latter alternative is also supported by the 
excellent agreement between model calculations and simulations 
(Field & Blackman 2002; Blackman & Brandenburg 2002). Fur- 
thermore, the simulations of Brandenburg & Subramanian (2005b) 
demonstrate that this quenching is accompanied by an increase of 
-«m toward ok, and that this «k itself is unquenched. Subsequent 
analysis of their data shows that ok remains unquenched regard- 
less of whether or not one uses the proper anisotropic expression 
(Brandenburg & Subramanian, unpublished). 

Our present results are of course restricted to the case of 
small magnetic and fluid Reynolds numbers. This means that we 
have not tested any of the actual closure assumptions, like the r- 
approximation. Such tests have so far only been done numerically 
(Brandenburg et al. 2004; Brandenburg & Subramanian 2005a). 
Clearly, at large magnetic and fluid Reynolds numbers the r- 
approximation can no longer yield exact results. Nevertheless, it 
provides a very practical tool to estimate the mean field transport 
coefficients in a way that captures correctly some of the effects that 
enter in the case of large magnetic and fluid Reynolds numbers. In 
that respect, it has been quite successful in reproducing the catas- 
trophic quenching result for periodic domains, as well as suggest- 
ing ways to alleviate such quenching. 

In Section 2 we outlined the conditions under which double 
FOSA is valid as being basically the requirement that R m <C 1 
and Re <JC 1. A subtle point concerns the validity of retaining the 
linear Lorentz force term, Bo Vb, while neglecting the nonlinear 
advection, u ■ Vtt, even though nonlinear advection is small com- 
pared to the viscous dissipation for Re <JC 1. This assumption is 
valid provided Bob/l > u 2 /I, or, using b ~ B m Bo, B 2 Rm > u 2 . 
(We thank Eric Blackman for pointing this out to us.) Note that in 
terms of the critical field B CT — ,/vfjqo, which divides the regimes 
where the Lorentz force is important (B > B cr ) and where it 
is not (Bo < Bcr), this requirement becomes Bo/B cr > Re 1 / 2 . 
Therefore, for small fluid Reynolds number our assumption of re- 
taining the linear Lorentz force term while dropping nonlinear ad- 
vection is indeed valid for most regimes of interest for a sup- 
pression. For smaller mean fields, where Bo/B cr < Re 1,/2 , in 
any case the Lorentz force has no impact. The interesting point 
seems to be that for low Reynolds number systems, the typical ref- 
erence mean field is not the equipartition field Bo = u, but rather 
B = B cr ~ u/(R m Re)^ 2 > u. 

Throughout this work we have adopted an externally imposed 
body force to drive the flow. This is commonly done in many sim- 
ulations in order to achieve homogeneous isotropic conditions that 
are amenable to analytic treatment. Clearly, this is not the case 
for many astrophysical flows that are driven by convection (e.g. in 
stars) or the magneto-rotational instability (e.g. in accretion discs). 
Such flows tend to show long-range spatial correlations, which 



© 2006 RAS, MNRAS 000, 1-13 



Kinetic and magnetic alpha effects 13 



means that the alpha tensor should really be treated as a integral 
kernel (see, e.g., Brandenburg & Sokoloff 2002). It is at present 
unclear whether such more natural forcings are closer to steady or 
to random forcing, and how big is the resulting (/, b) correlation. 
Given that this correlation represents already a linear effect, it is 
likely that the ap term can simply be subsumed into an expression 
for a modified kinetic ax. If this is the case, we can continue to 
write a « q?k + »m as the sum of a mildly suppressed kinetic part 
depending on the velocity field and a magnetic part, so that their 
sum accounts for the tendency toward catastrophic a-quenching in 
the absence of helicity fluxes. We recall that such a split is exact in 
the limit of delta-correlated forcing. 



6 CONCLUSIONS 

Our work was motivated in part by the detailed criticism expressed 
by Radler & Rheinhardt (2007). In view of our new results we 
can now make the following statements for small magnetic and 
fluid Reynolds numbers. Firstly, it is true that the qk that is calcu- 
lated under FOSA does indeed capture the full nonlinear a-effect 
— provided it is based on the actual velocity field. Secondly, the 
qk that is calculated in the r approximation, is not simply re- 
lated to «k, except in the limit of steady forcing. Thirdly, in the 
limit of small magnetic and fluid Reynolds numbers, both FOSA 
and the r-approximation give identical results. Indeed, all meth- 
ods of calculating the a-effect agree, as they should, given that the 
starting equations were the same. However, the force-field corre- 
lation cannot be ignored in general. The exception is when one 
analyzes the case of delta-correlated forcing in a manner akin to 
the r-approximation, where the force-field correlation does vanish 
and hence ckf = explicitly. In this case one can indeed write 
a — qk + am, or the sum of a kinetic and magnetic alpha ef- 
fects. Furthermore, due to the spatial non-locality of the Greens 
function for small magnetic and fluid Reynolds numbers, the r- 
approximation should be carried out at the level of spectral corre- 
lation tensors, as is done here. Somewhat surprisingly, the delta- 
correlated forcing case yields an asymptotic a oc B^ 2 scaling as 
opposed to the well-known a oc B^ 3 behaviour for steady forcing. 

Although our work is limited to small magnetic and fluid 
Reynolds numbers, the calculations of Method C and Method D, 
and its agreement with the results of Method A, (for both steady and 
time-dependent forcings), do suggest one way of thinking about the 
effect of Lorentz forces: they lead to a decrease of a predominantly 
by addition of terms proportional to the helical parts of the small 
scale magnetic field. Hence getting rid of such small scale mag- 
netic helicity by corresponding helicity fluxes, may indeed be the 
way astrophysical dynamos avoid catastrophic quenching of a to 
make their dynamos work efficiently. 
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